Probing CP violation with the electric dipole moment of atomic mercury 
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The electric dipole moment of atomic Hg induced by the nuclear Schiff moment and tensor- 
pseudotensor electron-nucleus interactions has been calculated. For this, we have developed and 
employed a novel method based on the relativistic coupled-cluster theory. The results of our theo- 
retical calculations combined with the latest experimental result of ^®^Hg electric dipole moment, 
provide new bounds on the T reversal or CP violation parameters ^qcd, the tensor- pseudotensor 
coupling constant Ct and (d„ — dd). This is the most accurate calculation of these parameters to 
date. We highlight the the crucial role of electron correlation effects in their interplay with the P,T 
violating interactions. Our results demonstrate substantial changes in the results of earlier calcu- 
lations of these parameters which can be attributed to the more accurate inclusion of important 
correlation effects in the present work. 
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The existence of a nonzero permanent electric dipole 
moment of a non-degenerate physical system is a sig- 
nature of the simultaneous violations of parity (P) and 
time-reversal (T) symmetries [T]. T violation implies the 
combined charge conjugation (C) and P violation ; i.e., 
CP violation via the CPT theorem 2]. The origin of 
CP violation is still not well understood. It has so far 
been observed only in the K j3|] and B [1, Q mesons and 
the results are essentially in agreement with the predic- 
tions of Kobayashi and Maskawa in the framework of the 
Standard Model 6]. However, this rnodcl cannot explain 
the matter-antimatter asymmetry [7| in the universe for 
which CP violation is a pre-requisite Q. In addition, it 
predicts atomic EDMs several orders of magnitude be- 
low their current limits Indeed, atomic EDMs 
are excellent probes of physics beyond the Standard 
Model d, and they provide important insights into 
a rich variety of CP violations-leptonic, semi-leptonic 
and hadronic sectors. Experimental searches are under- 
way for the EDMs of para-magnetic (open-shell) and dia- 
magnetic (closed-shell) atoms [ll|, ■ The results of the 
experiments can be combined with those of sophisticated 
atomic many-body calculations to determine various CP 
violating coupling constants at the atomic level which 
can ultimately be related to the CP violation parame- 
ters at the elementary particle level [l^. The EDM of 
dia-magnetic atoms arises predominantly from the nu- 
clear Schiff moment (NSM) and/or the electron-nuclcon 
tensor-pseudotensor interactions . These in turn arise 
from the nucleon-nucleon interactions or the EDM of nu- 
cleons, which originate due to the quark-quark interac- 
tions, EDMs and chromo EDMs of quarks at the elemen- 
tary particle level. 

In the present work, we concentrate only on the EDM 
of mercury (-'^^^Hg), a closed-shell atom. The present lim- 
its on important CP violation parameters like, Oqcd for 
strong interactions and the chromo EDMs of quarks have 
been obtained from the EDM of ^^^Hg pjj, which is the 



most accurate of all the data from atomic EDMs to date. 
The focus of our work in this Letter is to improve the cur- 
rent limits on the coupling constants associated with the 
electron-nucleon tensor-pseudotensor (T-PT) interaction 
(Ct) and the NSM(S'). A nonzero value of Ct implies 
physics beyond the Standard Model. The dependence 
of the T-PT interactions and the NSM on the nuclear 
spin makes closed-shell atoms, in particular, those hav- 
ing nonzero nuclear spin the best candidates to measure 
EDMs sensitive to the nuclear sector. 

For heavy atoms like mercury, it is customary to use 
the Dirac-Coulomb Hamiltonian i/ocj in atomic units 



Hdc = X! t'^"* ' + Awc^ + VN{ri)] + ^ 
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where refers to the electron coordinates, a and /3, 
the Dirac matrices and Vjv(''i), the nuclear potential, 
and the last term is the electron-electron Coulomb in- 
teraction. The key and the most challenging step in 
atomic many-body physics is to incorporate the effects 
of electron-electron Coulomb interaction, last term in 
Eq.([T]), as accurately as possible. Under the independent 
particle and the central field approximations [l4!|, with 
the introduction of the Dirac-Fock potential, the above 
Hamiltonian can be separated into an exactly solvable 
part (Hq) and a residual interaction part which consists 
of the Coulomb interaction and the Dirac-Fock potential 
[l^ . The single particle wave functions are computed 
self-consistently from Hq and the many particle wave- 
functions are expressed as Slater determinants built out 
of the single particle wavefunctions. The effects of the 
residual Coulomb interaction are calculated with many- 
body methods. For this, finite order many-body per- 
turbation theory (MBPT) and the configuration inter- 
action (CI) approach are two widely used methods p^ . 
An even superior method, with strong theoretical many- 
body physics underpinnings is the coupled-cluster theory. 
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In the present work, the P and T violating interactions 
are treated perturbatively to first order in addition to the 
electron-electron residual Coulomb interaction within the 
framework of the relativistic coupled-cluster theory. 

The wavefunction in coupled-cluster theory has an ex- 
ponential structure, see [l5| for a recent review of the 
method, to describe correlation effects in many-body sys- 
tems. It is non-perturbative and defines atomic states as 
superposition of states of different levels of excitations, 
which arise from the residual Coulomb interaction, with 
respect to a reference state. 

Mathematically, the coupled-cluster wavefunction can 
be expressed as 



I*, 



(2) 



where is the reference state containing a fixed num- 
ber of electrons and T^"^ is an operator which excites 
electrons out of it, thereby giving rise to states with 
different levels of excitations corresponding to differ- 
ent many-body effects. In our calculations we use the 
coupled-cluster singles and doubles approximation, that 
is r'^°) = T^^^ + T^^ . In second quantized form 

= ^ata,iP(0) and T^^ = ^ alala^ajZi^) , 



a.b,p,q 



excite one and two electrons respectively from the refer- 
ence state. The equations that determine the amplitudes 
of r*^°) are a set of coupled non-linear algebraic equations 
and these are solved iteratively till convergence. 

For closed-shell atoms, as mentioned earlier, one 
prominent source of EDMs is the nuclear Schiff moment 
S (NSM) , a P and T odd electromagnetic moment of the 
nucleus. For a finite size nucleus of radius Rn, the Schiff 
moment potential [16] is 
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where B = J p{R)R^dR and p{R) is the nuclear density. 
This potential interacts electrostatically with the elec- 
trons, it mixes atomic states of opposite parities and gen- 
erates a finite atomic EDM, cIa- Then, the atomic Hamil- 
tonian is i?atom = -ffoc+A^fpTV, where iJl^jiv^ = 
is the P and T violating interaction Hamiltonian and A is 
a T or CP violation parameter which can be considered 
as perturbation parameter. The eigenstates of the iJatom 
are the mixed parity states |^). To incorporate iJpTV 
as a first order perturbation, the exponential operator in 
coupled-cluster theory is redefined as e"^' ^+^'^^ The 
cluster operator T^^^ has one order of iJpTV and mixes 
the states of opposite parities. As a result of this, the 
ground state 
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Then, unlike the T^^^ equations, since -ffpTV is considered 
to first order only, the equations for the amplitudes of 



are a set of linear algebraic equations. 
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O — e^'' ' Oe^' ' where O is a general operator, is the 
normal-ordered Hamiltonian and |<I>q) are opposite parity 
Slater determinants. Further, like in the unperturbed 
cluster operators T^^\ we use the approximation T^^-* — 



T^^'' + T^"' . Then, the atomic EDM of the ground state 
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where D is the electric dipole operator. In the above 
expression, after expanding in terms of the cluster op- 
erators T(°) and T(^\ only the terms first order in T^^^ 
contribute. Often, c?a is computed perturbatively with 
the sum over states approach, which necessitates a trun- 
cation after the first few intermediate states. On the 
contrary, our relativistic coupled-cluster scheme doesn't 
involve summing over states and subsumes all possible in- 
termediate states within the chosen configuration space. 

Besides the NSM, the other possible source of EDM 
in closed shell atoms is the tensor-pseudotensor electron- 
nucleus interaction 



T-PT 
PTV 
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where Gf is fermi constant, Ct is coupling constant, 
(Ttv is nuclear spin and 7^ is Dirac matrix. It must be 
emphasized that, this form of interaction does not exist 
within the Standard Model of particle physics and Ct is 
zero. However, there are models which predict such an 
interaction (l3| . 

To extract the T or CP violation parameters, the 
atomic theory calculations are combined with the exper- 
imental data. In this context it is appropriate to rewrite 
Eq.® as 



(8) 



where r/ is the atomic enhancement factor. As defined 
earlier, the constant A is a T or CP violation parameter 
considered as a perturbation parameter. It can for ex- 
ample be the nuclear Schiff moment S or the coupling 
constant Ct- A precision atomic many-body calculation, 
like the coupled-cluster calculation reported here, would 
provide the value for a particular 77. Experimentally, the 
measured atomic EDM c?a is the sum total of contribu- 
tions from all the P and T symmetry violating phenom- 
ena within the atom. A bound on A = oJa/?? is obtained 
by combining the results of atomic theory and experi- 
mental data. Depending on the choice of the atom, it 
is possible to identify the dominant sources of T or CP 
violation and derive tighter bounds. 

For the present set of calculations, we ernploy the even- 
tempered Gaussian basis set expansion [13, [iBl- The 
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orbital basis set consist of (1 — 18)s, (2 — 18)pi/2,3/2, 
(3 - 13)^3/2,5/2, (4 - ll)/5/2,7/2, (5 - 8)37/2,9/2 and 
(6 — 7)/i9/2,ii/2- This orbital basis set is considered com- 
plete for the coupled perturbed Hartree-Fock (CPHF) 
calculations. That is, further increase in the number 
of orbitals does not change the results. In addition, we 
compute the ground state dipole scalar polarizability for 
i99jjg -y^g obtain a value 33.294eaQ, where oq is the Bohr 
radius, which is in excellent agreement with its experi- 
mental value [H, [10, im . 

To date, among the closed-shell atoms, ^^^Hg as men- 
tioned earlier sets the standard for the most precise EDM 
results. In a recent paper the new upper limit is re- 
ported as 

|d(i^3Hg)| < 3.1 X lO-'^e cm (95% C. L.). (9) 

Our atomic calculation based on the relativistic coupled- 
cluster theory gives 



jSchiff /199 
"■A 



(i99Hg) = -5.07 X 10" 



fm'^ 



e cm. 
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This is the first ever relativistic coupled-cluster result 
for any atomic EDM calculation arising from the NSM. 
Combining with the experimental result, the limit on the 
NSM is 



S'(^'^^Hg) < 6.1 X lO^^'^e fm 



(11) 



There is a large change of 96% from the result of coupled 
perturbed Hartree-Fock calculation 
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Hg) = 2.8 X 10 
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fm3 



e cm. 
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reported earlier |22l | , which is in excellent agreement with 
the result of a similar calculation but in the framework 
of the relativistic coupled-cluster theory [2^ . The large 
change in the two results demonstrates the importance 
of electron correlation effects and their interplay with the 
Hp'^y^ interaction in determining the magnitude of the 
NSM. 

It is possible to separate the contributions of individ- 
ual terms in Eq.® and the many-body perturbation di- 
agrammatic representation of the dominant terms are 
shown in Fig[T] These diagrams represent the excita- 
tions and de-excitations due to cluster operators and 
dressed dipole operator D. Earlier calculations [H, H^l 
incorporate only a certain class of two-particle two-hole 
excitations which are subset of the correlation effects 
we have included through the cluster operator T^^^ in 
the present calculation. The two most dominant terms 



are T^^^^ D and T^'''DT^"', in Eq.dTO]) these terms in- 
dividually contribute —5.40 x 10~^^(5'/(e fm^))e cm and 
-0.17 X 10~^'^{S/{e fm3))e cm respectively. In FigUl the 
diagrammatic equivalent of these terms are (a) and (b) 
respectively. 

Our result of ^^^Hg EDM arising from the electron- 
nucleus tensor-pseudotensor interaction is 



.(i)tr,^(0) 



tT-PT 



4.3 X lO^^^CrCTNe cm. 



(13) 
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FIG. 1: Two of the dominant many-body perturbation dia- 
grams of dA- Diagrams (a) and (b) represent the terms Tj^'^D 

and T^'^^^DTj"'' respectively. At the interaction vertices, the 
horizontal wavy blue, black and red lines represent the HpTv, 
Tj"' and the electric dipole operator D respectively. 



Commred to the CPHF result -6.19 x IQ-'^^CrcrNe cm 
[HiH^I, the change with the additional correlation effects 
is not so dramatic. There is a decrease of 31%, which 
is significant but not so spectacular as in d^^^^ {^^^Hg). 
This comparison demonstrates, without any ambiguity, 
the importance of electron correlation effects in precision 
atomic EDM calculations. A closer examination of the 
structure of the two P and T violating Hamiltonians in 
the present work sheds some light on why the electron 
correlation effects are larger in the case of the NSM than 
the tensor-pseudotensor interaction. Both the pi^2 and 
P3/2 electrons are actively involved in the interplay of the 
i/pTV and electron correlation effects in the former, while 
the contribution of the P3/2 electrons is negligible in the 
latter [13]. 

The individual contributions follow similar trend as 
in the case of NSM. The terms T^^'>''d and T'-^^''dT^°^ 
give the largest (« 95%) and the second largest contri- 
butions to d 



rT-PT 

A ' 



4.8 X 10-^^CT<JNe cm and -0.27 x 



10 Cto'nb cm respectively. Then, a limit 
Ct < 1.4 X 10"^ 



(14) 



is obtained after combining our results with the experi- 
mental data. 

Assuming that the NSM arises from the nucleon- 
nucleon interactions with pions as the dominant medi- 
ators [li 



- 0.075!'^^ 



+0.025^^ 



i(2) 



NN 



■ fm^ 



(15) 

where g-„NN and 5^*^^ are the CP conserving and CP 
violating pion-nucleon coupling constants respectively. 
Here, i — 0,1, and 2 represent isoscalar, isovector and 
isotensor respectively. Considering g^l;^^ as the most 
dominant 



< 6.4 X 10-13. 



(16) 
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The coupling constant <?jr_/VAr related to the chromo- 
EDMs of quarks [1^ , from the above result 

{du - dd) < 3.2 X IQ-'^'^e cm. (17) 

Next, consider the maximum contribution to NSM aris- 
ing from g^"lj^, then 

5i°^jv < 4.5 X 10-12. (18) 

Since, g^^^j^ = 0.0276Iqcd 113, we get the bound 

0QCD < 1-7 X 10-10. (19) 

The value we have obtained for the NSM is likely to 
give the most stringent bounds for supersymmetric CP 
violating phases [l^, [HI. In addition, from our re- 
sults and the experimental data, it is also possible to set 
improved limits on the specific CP violating parameters 
predicted by various extensions of the Standard Model, 

gSUSY^ eHiggs^^LR ^ 

In conclusion, we have developed a unique relativistic 
coupled-cluster based many-body method that takes into 



account the physical effects arising from the interplay of 
two very different kinds of fundamental interactions-the 
CP conserving electron-electron Coulomb and CP violat- 
ing electron- nucleus interactions. The results obtained 
for the EDM of ^^^Hg by the application of this method 
and the latest experiment on this atom fll] yield the 
most accurate limits to date on some important CP vio- 
lating parameters. The electron correlation effects play a 
critical role in improving the existing limit on these pa- 
rameters. These limits constrain the possible extensions 
to the Standard Model, thereby enhancing our current 
knowledge of the intriguing phenomenon of CP violation. 
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